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DUCTS  FOR  ACCELERATED  FLOW1 
by 

T.  C.  Lin2 


ABSTRACT 

A method  of  design  for  two-dimensional  contracting 
ducts  for  incompressible  and  inviscid  flow  has  been  developed 
which  gives  a rapid  change  of  cross-section  with  no  pressure 
discontinuity  or  adverse  pressure  gradient.  While  giving 
essentially  uniform  flow  on  both  ends  it  yields  ducts  of  smaller 
overall  length  and  shorter  wall  length  than  those  given  by 
existing  methods  for  any  contraction  ratio.  The  design  is 
based  on  a fan- shaped  boundary  chosen  in  the  hodograph  plane, 
for  which  the  flow  is  found  by  conformal  mapping.  The  ducts 
have  a longer  inlet  portion,  part  of  which  may  be  used  as  a 
settling  chamber,  a short  central  portion  for  the  acceleration 
of  the  flow  along  the  wall,  and  a short  outlet  portion,  thus 
cutting  down  power  losses.  As  examples  of  the  method  the 
design  of  ducts  of  contraction  ratios  3 and  10  is  carried  out. 


1.  The  results  presented  in  this  paper  were  obtained  in  the 

course  of  research  sponsored  by  the  Office  of  Naval  Research 
under  Contract  N7onr-35801  with  Brown  University, 

• Research  Associate,  Division  of  Applied  Mathematics,  Brown 
University,  Providence  12,  R,  I, 
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l.  Introduction  and  Summary. 

It  is  of  interest  in  the  field  of  fluid  mechanics  to 
determine  the  shape  of  the  contracting  duct  that  will  most  effec- 
tively increase  the  velocity  of  a stream  of  fluid  by  a ratio  r. 

In  order  to  reduce  the  boundary  layer  effect  the  wall  length  of 
a contracting  duct  should  be  as  short  as  possible,  while  keeping 
essentially  uniform  flow  on  both  ends  of  the  duct.  In  addition 
these  should  be  no  pressure  discontinuity  or  adverse  pressure 
gradient  along  the  wall  to  cause  separation  and  unsteady  conditions 

In  the  present  investigation  two-dimensional,  irrotational 
flow  is  considered.  The  fluid  is  assumed  to  be  incompressible  and 
inviscid.  Hor  a real  fluid  the  shapes  of  the  contracting  ducts 
should  be  modified  by  a slight  displacement  of  the  wall  outward 
to  allow  for  the  boundary  layer  thickness.  In  the  case  of  a com- 
pressible medium  the  compressibility  effect  on  the  shapes  can  be 
obtained  by  using  the  Karman-Tsien  linearized  pres sure -volume 
relation. 

Methods  of  design  for  contracting  ducts  in  two-dimension^ 
irrotational,  incompressible  and  non-viscous  flow  have  been  given 
by  various  investigators,  notably,  Hughes  [1]*,  Cheers  [2], 
Lighthill  [3],  Whitehead,  Wu  and  Waters  [4],  and  Libbey  and  Riess 

m. 

The  method  by  Hughes  uses  the  isobaric  lines  as  the  wall 

of  the  duct  and  obtains  the  most  rapid  change  of  cross-section 

possible  without  reversal  of  pressure  gradient,  but  it  introduces 

a spiral  singularity  and  hence  a pressure  discontinuity  along  the 

* Numbers  in  square  brackets  refer  to  Bibliography  at  end  of 
paper. 
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wall,  and  has  to  be  modified  suitably  for  practical  applications. 
The  others  avoid  the  discontinuity  defect  but  fail  to  make  use  of 
the  isobaric  lines  and  hence  give  longer  ducts. 

Besides,  all  these  methods  have  the  drawback  that  the 
maximum  angle  of  the  wall  of  the  duct  is  not  adjustable.  In 
Hughes'  case  the  angle  of  the  wall  of  the  duct  tends  to  infinity 
at  the  singular  point  for  all  contraction  ratios  and  causes  the 
pressure  discontinuity,  while  in  all  the  other  cases  for  a given 
contraction  ratio  there  is  a fixed  maximum  wall  angle,  which  is 
either  too  small  or  too  large  to  give  an  optimum  shape  and  leads, 
to  ducts  longer  than  necessary. 

In  the  present  paper  a method  of  two-dimensional  design 
for  contracting  ducts  is  given.  Apart  from  fully  utilizing  the 
isobaric  lines  for  the  inlet  and  outlet  portions  of  the  wall,  it 
introduces  a parameter  a for  adjusting  the  maximum  wall  angle  and 
eliminating  the  spiral  singularity.  The  advantage  of  using  the 
isobaric  lines  is  explained  in  section  2.  For  a given  contraction 
ratio  the  optimum  a is  obtained  by  minimizing  either  the  overall 
length  of  the  duct  or  the  length  of  the  wall  of  the  duct.  In 
order  to  reduce  the  thickness  of  the  boundary  layer,  ducts  of 
the  shortest  wall  length  are  of  primary  consideration. 

While  giving  essentially  uniform  flow  on  both  ends  and 
no  pressure  discontinuity  or  adverse  pressure  gradient  along  the 
wall,  the  present  method  yields  ducts  of  smaller  overall  length 
and  shortest  wall  length  than  those  given  by  existing  methods  for 
any  contraction  ratio.  Besides,  the  ducts  have  a longer  inlet 
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portion,  part  of  which  may  be  used  as  a settling  chamber,  a short 
central  portion  for  the  acceleration  of  the  flow  along  the  wall, 
and  a short  outlet  portion,  thus  cutting  down  power  losses. 
Typical  designs  are  presented  for  ducts  of  contraction  ratios  3 
and  10  as  examples  of  the  methods. 


Notations. 

c,c’ 

fl  '*•  f9 
F(k,cp) 

i 


constants,  Eqs,  (25)  and  (32) 

functions,  Eqs.  (11)  to  (14),  (20)  to  (22),  and  (40) 
Legendre  elliptic  integral  of  the  first  kind,  Eq.  (4) 
t/Ti 


Im 

k 

k» 

K 

K» 

r 

Re 

s 

S 

t 

u,v 

U 

w 

x,y 

z 

a 

«o 

Y 


Imaginary  part  of 
parameter,  Eq.  (2) 
parameter,  Eq.  (4) 

F(k,  m/2),  Eq.  (4) 

F(lc»,  x/2),  Eq.  (4) 
contraction  ratio,  Eq.  (3) 
real  part  of 

length  of  the  wall  of  the  duct,  Eq.  (37) 

iQ 

= pe  , complex  variable,  Eq,  (2) 
complex  variable,  Eq.  (1) 

velocity  components  in  the  x and  y-directions 

velocity  upstream  of  the  duct,  Eq,  (1) 

= u 7 *v>  complex  variable  in  the  hodograph  plane, 

Eq.  (1)  ’ 

Cartesian  co-ordinates  in  the  physical  plane,  Eq.  (6) 
= x + iy 

maximum  wall  angle,  Eq.  (2) 
optimum  a 

sin"1  k 
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r* 

42 

vl»  V2 
0 
P 

<P 

*1 

Q 


sin"1  k* 

velocity  deviations  at  the  upstream  and  downstream 
ends  of  the  duct,  Eqs.  (28)  and  (35) 

permissible  velocity  deviations  at  the  unstream 
and  downstream  ends  of  the  duct,  Eq.  (49) 

angular  coordinate  in  the  S-plane,  Eq.  (2) 

radial  distance  from  the  origin  in  the  S-plane, 

Eq.  (2) 

parameter 

tan“1(k1//2) , Eq.  (37) 

parameter,  Eq.  (41) 

complex  potential  function,  Eq.  (5) 


2.  The  Boundary  Shape  An  the  Hodograph  Plane. 

The  shape  of  the  boundary  of  the  duct  is  first  chosen 
in  the  hodograph  plane,  where  the  co-ordinates  are  the  velocity 
components  u and  -v.  This  is  most  convenient,  for  in  the  hodo- 
graph plane,  concentric  circles  centered  at  the  origin  are  isobaric 
or  constant  velocity  lines,  and  straight  lines  through  the  origin 
are  lines  of  constant  flow  direction. 

Assume  that  the  duct  is  symmetrical  with  respect  to  the 
x-axis,  Fig.  1,  and  has  an  upstream  velocity  U,  and  that  the  con- 
traction ratio  is  r.  In  the  hodograph  plane,  Fig.  2,  the  up- 
stream end  of  the  duct  corresponds  to  the  point  A,  (U,0),  and  the 
downstream  end  to  the  point  F,  (rU,0),  All  the  stream  lines 
start  from  A and  end  in  F.  The  axis  of  symmetry  in  the  physical 
plane  corresponds  to  the  straight  line  AF  in  the  hodograph  plane. 

In  order  to  avoid  an  adverse  pressure  gradient  along 
the  wall  of  the  duct,  the  boundary  in  the  hodograph  plane  must 
lie  on  or  between  the  two  concentric  circles  of  radii  U and  rU 
centered  at  the  origin  0,  If  a is  the  maximum  angle  between  the 
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wall  and  the  center  line  of  the  duct,  the  boundary  in  the  hodo- 
graph  plane  must  lie  on  or  between  the  two  straight  lines  passing 
through  the  origin  0 and  making  angles  of  +a  with  the  u-axis. 
Therefore  the  boundary  in  the  hodograph  plane  must  lie  within 
the  fan-shaped  domain  or  on  its  boundary.  The  upper  half  of  the 
fan-shaped  domain  is  shown  in  Fig,  2, 

From  the  symmetry  of  the  problem  it  is  only  necessary 
to  consider  the  upper  half  of  the  duct.  The  wall  of  the  duct  is 
conveniently  divided  into  three  portions,  the  inlet  BC,  the 
central  portion  CD,  and  the  outlet  DE,  Fig.  1,  For  a short  inlet 
BC  it  is  desirable  to  turn  the  flow  along  the  inlet  portion  of  the 
wall  to  the  maximum  wall  angle  a as  fast  as  possible  without  in- 
creasing the  speed  of  the  flow.  For  a rapid  contraction  at  the 
central  portion  it  is  desirable  to  increase  the  speed  of  the  flow 
along  the  wall  at  this  maximum  angle  a as  much  as  possible  up  to 
the  final  speed  rU,  For  a short  outlet  DE  it  is  desirable  to  turn 
the  flow  along  the  outlet  portion  of  the  wall  back  to  the  horizontal 
direction  as  fast  as  possible.  It  follows  from  these  physical 
arguments  that  for  ducts  of  smaller  wall  length  the  boundary  in 
the  hodograph  plane  should  coincide  with  the  boundary  of  the  fan- 
shaped domain,  Fig.  2. 

As  a comparison  the  boundaries  used  by  Whitehead,  Wu 
and  Waters,  and  by  Libby  and  Riess  are  also  shown  in  Fig.  2,  The 
extreme  case  that  a approaches  infinity  gives  the  boundary  used 
by  Hughes.  The  boundaries  used  by  Cheers  and  Lighthill  are  com- 
pared with  the  present  boundary  in  Fig.  3. 
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3.  The  Mapping  Functions. 

The  flow  in  the  hodograph  plane  is  found  by  conformal 
mapping.  Let  w = u - iv  be  the  complex  co-ordinate  in  the  hodo- 
graph plane.  The  upper  half  of  the  fan-shaped  domain  in  the 
w-plane,  Fig,  2,  is  mapped  into  the  rectangule  in  the  t-plane, 
Fig.  3j  by  the  equation 

w = U et1  (1) 

where  1 < |w|  AJ  < r,  0 < Arg  w < a}  and  0 < Re  t < In  r, 

0 ^ Im  t < a. 

The  rectangle  in  the  t-plane  is  mapped  into  the  upper 
half  of  the  S-plane,  Fig.  4,  by  the  equation, 


t = JL  [S  as 

2K  Jo  /sTs  +'  k) (kS  + 1) 


(2) 


where  Im  S > 0.  The  parameter  k is  related  to  the  contraction 
ratio  r and  the  maximum  wall  angle  a by  the  equation 


In  r = ctK’A  (3) 

where  K and  K*  are  complete  Legendre  elliptic  integrals  of  the 
first  kind  defined  by 

K = F(k,ii/2), 

K»  = F ( k ’ , n/2 ) , 

*>„).[»■— ffl. . ( (4) 

J 0 ^1~-  k^sin^cp 

k»2  = 1 - k2. 

It  is  readily  seen  that  K*/K  is  a function  of  k or  k1  only.  The 
values  of  K»/K  for  various  values  of  k or  k1  are  given  in  tables 
of  complete  elliptic  integrals. 
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The  complex  potential  function  for  the  upper  half  of 
the  S-plane  with  a source  at  the  origin  and  a sink  at  infinity  is 
simply 

Q = (U/n)  In  S.  (5) 


The  strength  of  the  source,  U A,  is  chosen  to  yield  a half  width 
of  unity  for  the  entrance  of  the  duct. 

If  z = x + iy  is  the  complex  co-ordinate  in  the  physical 
plane,  one  has 

z-fsS 

vj  w 

Eliminating  w,  t,  and  Q between  Eqs.  (1),  (2),  (5),  and 
z = - ( Exp  [ - JL  \ S dS ] , 

J 1 A s 2K  J o /s  ( S + k)  fkS  + 1) 


(6) 

(6)  yields 
(7) 


The  lower  limit  lA  of  the  above  integral  is  chosen  for  mathe- 
matical simplicity  and  to  fix  the  origin  of  the  z-plane  on  the 
center  line  and  close  to  the  throat  of  the  duct. 


The  Boundary  Shane  in  the  Physical  Plane. 

As  S varies  from  zero  to  minus  infinity  along  the  nega- 
tive real  axis,  z traces  along  one  branch  of  the  boundary  of  the 
duct  from  the  upstream  end  to  the  downstream  end.  Putting  S = - p, 
using  the  appropriate  limits,  and  separating  the  real  and  imaginary 
parts,  one  obtains  from  Eq.  (7)  the  following  equations  for  the 
boundary  of  the  duct  in  the  physical  plane  in  terms  of  the 
parameter  9 . 

For  the  inlet  BCs 

x - 3^  = -[f^Ckja,*)  - f2(k,  a,  SP)3, 
y - yc  = f3(k,a,<p),  j 


\ 
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where  0 < 9 < n/2,  and  cos2  cp  = p/k. 

For  the  central  straight  portion  CD: 

■s 

Xq  - XD  = - f^Ck'  ,r,  n/2)cos  cc, 
yc  - yD  = f^Ck*  ,r,Tt/2)sin  a. 


For  the  outlet  DE: 


- x = - J [f'1(k,  a,  cp 

- y = J (k,  a,cp) , 


a, 9)  - f2(k,a,  9)]  , 


(10) 


where  0 < 9 < tc/2,  and  sec  9 = kp. 

A comparison  of  Eqs,  (8)  and  (10)  shows  that  the  inlet 
and  outlet  portions  of  the  wall  follow  the  same  geometrical  curve 
and  differ  only  by  a scale  factor  1/r.  The  half  widths  of  the 
duct  at  the  upstream  end  and  the  downstream  end  approach  asymp- 
totically to  1 and  1/r,  respectively. 

In  the  above  the  functions  f's  are  defined  as  follows: 

f-,(k,a,  9)  = 2 cos[aF(k,  £ - 9 )/K]ln  sec  9,  (11) 

2a  f 9 sin  [aF(k,  5-9  )/K]ln  sec  9 d9 

f2(k,a,<p)  = — | %—  -- _■  ■ T72 * (12) 

,|  0 (1  - k2  cos2?)1' 


f (k,a,9)  = - sin[aF(k,  ^ - 9)/K]tan9d9  , (13) 

0 71  JO 

of  9 -[1  - F(k«,  § -9)/K»]  . _ . 

_ „ , \ u2  Y r ’ 2 sin  29  d9 

V1 *’>'.*>  ■ V L r — i 1 — OH) 

^0  1 - k*2  cos^  9 


These  integrals  can  be  evaluated  without  difficulty. 
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5.  The  Velocity  Distribution. 

Combining  Eqs,  (1)  and  (2)  gives  the  expression  for  the 
velocity  in  terms  of  S: 


w = U ExpCJL 
2K 


dS 


0 /S(S  + k)  (kS  + l) 


],  S = pe19, 


(15) 


where  Im  S > 0,  1 < Iwl/Q  AJ  £ r,  and  0 < Arg  w £ a. 

Consider  first  the  velocity  distribution  along  the  wall 
of  the  duct.  From  the  boundary  chosen  in  the  hodograph  plane  it 
is  clearly  seen  that  along  the  inlet  BC 

I w | = U,  (16) 


and  along  the  outlet  DE 


I w | = rU.  (17) 

Therefore  along  the  wall  of  the  duct  the  acceleration  of  the  fluid 
takes  place  only  along  the  central  straight  portion  CD,  Using 
the  condition  © = m,  one  obtains  from  Eqs.  (15)  and  (7)  the 
velocity  distribution  along  CD  in  terms  of  the  parameter  cp  as 
follows : 


(18) 


|w|/U  = f^(k* , r , cp) , 1 

| z - zc  | = f^  ( k * , r , 9 ) , J 

where  0^9^  m/2,  and  cos2  9=  (1  - kp)/(l  - k2). 

Similarly,  using  the  condition  © = 0,  one  obtains  the 
velocity  distribution  along  the  center  line  of  the  duct, 

u/U  = f ej(k'  ,r,  9) , 
x = f6(k»,r,9)  - f (k«,r,9), 

where  0 9 £ it/2,  and  tan2  9 = kp.  From  symmetry,  v is  zero 

along  the  center  line. 


(19) 
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In  the  above,  f^  Is  defined  by  Eq.  (14);  f^,  fg,  and  f^ 
are  defined  as  follows: 


f5(k.,r)9)  = rtl-F(k’-!-^,;1 


(20) 


f6(k«,r,cp)  = | r"^1  ‘ F(k'»  2 “ ^^ln  tan  <p , (21) 


f?(k*,r,cp)  = 


2 In  r 
ftK» 


cp  -[1-F(k»,  | - tp  )/K»  3 


tcA 


(1  - k»2  cos2  cp  )1^2 


In  tan  9 d(p  ^2) 


Consider  next  the  velocity  distribution  at  the  upstream 
end  of  the  duct.  At  the  upstream  end 

|S|/k  « 1.  (23) 

On  neglecting  sA  and  its  higher  powers  as  compared  to  1,  and 
since  o/K  is  of  the  order  1 or  less,  Eqs.  (15)  and  (7)  reduce, 
respectively,  to 

w/U  = 1 - (aA)(S/k)1/2,  (24) 

and 

z = - ~[ln(k/S)  + 2(aA)(sA)1/2]  + c.  (25) 

ft 

Here  c is  a real  constant  and  its  value,  in  general,  depends  on 
the  two  parameters  k and  a . 

Separating  the  real  and  imaginary  parts  of  Eqs,  (24)  and 
(25)  yields 

uAJ  = 1 - (a/K)(p  A)1//2  cos  (9/2), 
v/U  =-(aA)(p  A)1/2  sin  (©/2),  ^ 

-J 

x = - i[ln(k/p)  + 2(a/K)  (p  A)-1^2  cos  (0/2)  ] + c, 

ft 

y = i[e  - 2(aA)(pA)1/2  sin  (0/2)], 


and 


(26) 


J (27) 
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where  p/k  <<  1,  and  0 £ | © | £ u.  From  symmetry  © is  negative  for 
the  lower  half  of  the  duct. 

Denoting  the  velocity  deviation  at  the  upstream  end  of 
the  duct  by 

JJL 1 = | w - u|  Al,  (28) 

which  is  a measure  of  the  ratio  of  the  undesirable  velocity  devia- 
tion to  the  uniform  velocity  at  the  upstream  end,  one  has  from 
Eq.  (24) 

ix  = (cc/K)(p/k)1/2,  (29) 

where  pA  « 1,  From  Eqs.  (27)  it  is  readily  seen  that  for 
constant  x the  maximum  p occurs  at  y = © = 0,  It  follows  from 
Eq.  (29)  that  at  the  upstream  end  of  the  duct  for  constant  x 
the  maximum  occurs  at  the  center  line  of  the  duct. 

Similarly,  at  the  down  stream  end  of  the  duct 

l/(k Is  I ) « 1,  (30) 

and  Eqs.  (15)  and  (7)  reduce  to 

w/rU  = 1 - (a/K)(kS)"1/2,  (31) 

and 

z = ± [ ln(kS)  + 2(a/K)(kS)“1/2]  + c*,  (32) 

where  c'  is  a real  constant.  Separating  the  real  and  imaginary 
parts  yields 

l/2 

u/rU  = 1 - (a/K) (kp)  cos  (©/2), 
v/rU  = -(a/K)(kpr1/2  sin  (©/2). 
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and 


, 1 /o 

x = ~[ln(kp)  + 2(a/K)(kp)~  cos(0/2)]  + c 

y = ~ [©  ~ 2(a/K)(kp)”1/2  sin  (0/2)], 


K3^> 


where  1/k  p « 1 and  0 < I 9 I < x. 

Similarly,  denoting  the  velocity  deviation  at  the  down- 
stream end  of  the  duct  by 


p2  = |w  -rU|/rU , (35) 

one  has  from  (31) 

jx2  = (aA)(kp)"1/2,  (36) 


where  lAp  <<  1.  From  Eq.  (34)  it  is  readily  seen  that  for  con- 
stant x,  the  minimum  p occurs  at  y = 9 = 0.  It  follows  from  Eq. 
(36)  that  at  the  downstream  end  of  the  duct  for  constant  x the 
maximum  occurs  at  the  center  line  of  the  duct. 

The  velocity  distribution  across  the  central  section  of 
the  duct  is  well  illustrated  by  the  isobaric  and  equal-velocity 
potential  line  LM,  Fig,  1.  The  line  LM  corresponds  to  the  upper 
half  of  the  unit  circle  in  the  S-plane,  Fig.  4. 


P = 


The  point  L lies  on  the  center-line  of  the  duct.  Putting 
1,  one  obtains  from  Eq.  (19) 


xL  = f^k'jr,^)  - f7(k«,r,?1) 
92  = tan-1(k1//2) 


(37) 


Along  the  line  LM p = 1;  Eqs. 
respectively, 


(15)  and  (?)  yield, 


| w|/U  = r1/2 


(38) 
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and 


where 


8 


X - 

■ XL 

-1/2 

= r 

f8(k,a, 

!<P) 

y 

= r-l/2 

l<P) 

= 2 

f 9 

sin[ 

a 

F(2id^ 

I 9 

)3d9 

ft 

J 0 

(l+k)K 

1+k 

= 2 

f 9 

\ 

cost 

a 

,2k1/2 
F( , 

. 9 

)3d9 

X 

j 0 

(l+k)K 

1+k  1 

» * 

(39) 


r (4o) 


0 < 9 = 4 < 2 j and  9 = 2 at  the  point  M.  On  using  Landen’s 

C*  C.  £— 

Transformation  (6], 

sin  2<p 


tan  ij)  = 


k + cos  2<p 


F(2irr,ip>  * F(M)’ 

since  tJj  = u for  9 = rc  /2,  it  follows  that 

1/2 


(41) 


F( 


2k 


■•?)-! 


(l+k)K  1+k  * 2 
The  distance  CM  is  obtained  from  Eq.  (18)  by  putting 


(42) 


P=  1. 


ZC  " ZM  ^ = *4^*1  r»  9j) 
9^  = tan"1(k1//2) 


(43) 


6.  The  Length  of  the  Wall  and  the  Optimum  g. 

Let  s denote  the  length  of  the  wall  of  the  duct.  Since 
on  the  wall  |©  | = x,  Eqs,  (5)  and  (6)  yield 


-s  is- 


pw 

Along  the  inlet  BC  Jw  j = U,  and  Eq,  (44)  gives 


(44) 


\ 
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SBC  = 5 ln<M=B>,  0*5) 

where  0 £pg/k  < 1.  Similarly,  along  the  outlet  DE  |w  | = rU,  and 

SDE  * it  (W 

where  0 ^ l/kpE  £ 1.  From  Eqs.  (9), 

sCD  = f4(k»,r,it/2).  (^7) 

Therefore  the  total  length  of  the  wall  from  B to  E is 

SBE  = i“P-n(k/pB)  + i In(kpg)  ] + f^(k' ,r,x/2),  (48) 


So  far  a is  still  undetermined,  but  it  is  related  to 
the  contraction  ratio  r and  the  parameter  k or  k1  by  Eqs.  (3)  and 
(4).  For  a given  r one  may  assume  a value  for  a and  find  the 
corresponding  value  of  k'  from  Eqs.  (3)  and  (4).  In  general, 
however,  for  a given  contraction  ratio  r with  a given  set  of 
permissible  velocity  deviations  v-^  and  at  the  two  ends  of  the 
duct,  there  is  an  optimum  a and  hence  an  optimum  k*  that  will 
give  a shortest  wall  length  for  the  duct.  This  is  most  desirable 
as  it  will  minimize  the  boundary  layer  effect  to  the  flow. 

In  order  to  determine  the  optimum  a or  k'  it  is  neces- 
sary to  express  the  length  of  the  wall  in  terms  of  the  contraction 
ratio  r,  the  permissible  velocity  deviations  and  and  a or 
k'.  Let  H and  G denote  respectively  the  two  points  at  the  up- 
stream and  downstream  ends  and  on  the  center  line  of  the  duct, 

Fig,  1,  Since  the  maximum  velocity  deviations  on  both  ends  of 
the  duct  always  occur  at  the  center  line,  it  follows  from  Eqs. 

(29)  and  (36),  that 


\ 
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v1  = (a/K)(pH/k)1/2, 

v 2 = (a/K)(PGk)“1/2. 


(49) 


Using  the  condition  that  the  value  of  9 is  zero  at  the  center  line 
and  it  at  the  upper  wall,  one  has  from  Eqs.  (27)  and  (34), 

—V 


ln(k/pB)  = ln(k/p  ) + 2(a/K) (p jj/k)1^2 


and 

which, 

and 


ln(kpE)  = ln(kpQ)  + 2(a/K)(kpQ)‘1/2, 
in  view  of  Eqs.  (49)  and  (3),  become 

ln(k/pB)  = 2[v1  + ln(^0] 


(5D 


ln(kp  ) = 2 [v  + ln(iS-£)]. 
t v^K  * 

Substituting  these  into  Eq.  (48)  gives  the  length  of  the  wall  in 
terms  of  r,  v^,  and  k1  as  follows. 


+ 


+ 


- (1  + i)ln  (r^— — ) ] + f,(k*  ,r,Tc/2). 
r In  r 4 7 ’ 


(52) 


Equation  (52)  holds  for  v ^ and  small,  which  is  the  case  of 
practical  interest.  Both  In  K*  and  f^(k*  ,r , n/2)  increase  with  kf, 
and  for  a given  r there  is  an  optimum  k’  and  hence  an  optimum  a 
which  gives  the  shortest  wall.  Minimizing  Eq.  (52)  one  obtains 
the  optimum  k1,  and  hence  from  Eqs.  (3)  and  (4)  the  optimum  a. 

The  latter  is  platted  against  r in  Fig,  5.  As  a comparison  the 
a’s  given  by  the  other  methods  are  also  shown  in  the  same  figure. 
The  variation  of  sin"1  k against  the  compression  ratio  r for 


\ 
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various  values  of  a is  shown  in  Fig.  5a«  The  variations  of  k, 
k',  K,  Kf,  K/K*  and  K’/K  against  y and  y*  are  shown  in  Fig.  5b, 
where  y=  sin"1  k and  y'  = sin"1  k*. 

7.  Examples. 

Calculation  of  the  shape  of  the  boundary  has  been  made 
for  the  contraction  ratios  r = 3 and  10  using  the  corresponding 
optimum  a’s.  The  results  are  given  in  Tables  1,  la,  and  2,  2a, 
and  are  illustrated  in  Figures  6,  6a,  6b  and  7,  7a,  7b,  respec- 
tively. 

8.  Conclusion  and  Discussion. 

Using  the  constant  velocity  lines  as  the  inlet  and 
outlet  portions  of  the  wall  and  introducing  a parameter  a for 
adjusting  the  maximum  wall  angle,  the  present  method  yields, 
for  all  contraction  ratios,  contracting  ducts  of  shorter  wall 
length  and  smaller  overall  length  than  those  given  by  existing 
methods  with  no  pressure  discontinuity  or  adverse  pressure  gradient. 

The  ducts  have  a longer  inlet  portion,  part  of  which  may 
be  used  as  a settling  chamber,  a short  central  portion  for  the 
acceleration  of  the  flow  along  the  wall , and  a short  outlet 
portion,  thus  cutting  down  power  losses. 

Physical  arguments  and  numerical  results  indicate  the 
wall  of  the  duct  given  by  the  present  method  to  be  the  shortest 
possible  without  pressure  discontinuity  or  adverse  pressure  grad- 


\ 
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The  present  investigation  may  be  extended  to  ducts  with 
a certain  allowable  adverse  pressure  gradient  or  ducts  with 
parallel  walls  on  one  or  both  ends  by  shifting  the  locations  of 
the  source  and  sink  in  the  S-plane  along  the  positive  real  axis. 

For  a compressible  medium  the  effect  of  compressibility 
on  the  shape  of  the  ducts  may  be  obtained  by  using  the  Karman-Tsien 
linearized  pressure-volume  relation. 

Boundary  shpaes  which  are  satisfactory  for  two-dimension- 
al ducts  will,  in  general,  have  some  extra  margin  of  safety  for 
the  pressure  gradient  when  used  for  tunnels  of  circular  cross 
section,  in  spite  of  the  larger  contraction  ratio  [4],  Making 
use  of  the  close  similarity  between  the  shapes  of  the  stream  line.7 
for  the  two  flows,  this  can  readily  be  checked  for  any  particular 
case  when  the  two-dimensional  solution  is  known.  In  comparison 
with  the  ducts  for  axi-symmetric  flow  given  by  Whitehead,  Wu,  and 
Waters  [4],  it  is  expected  that  the  present  design  when  used  for 
tunnels  of  circular  cross-section  will  also  offer  a shorter  wall. 

Flow  through  Y-shaped  channels  with  the  same  contraction 
ratio  as  the  ducts  and  with  no  adverse  pressure  gradient  can  be 
obtained  by  shifting  the  source  from  the  point  A to  C in  the 
S-plane,  and  making  use  of  the  symmetry  of  the  problem.  One  half 
of  the  Y-shaped  channels  yields  the  elbows  for  accelerated  flow 
given  by  Carrier  [7]. 

Similarly,  flow  through  inverted  Y-shaped  channels  can 
be  obtained  by  keeping  the  source  at  A and  shifting  the  sink  from 
F to  D in  the  S-plane.  This  method  can  be  generalized  to  study 
flow  through  an  intersection  of  several  channels  of  various 
velocity  ratios  and  at  various  angles. 


\ 
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TABLE  1 

Values  of  f^,  f2,  f ^ , fg  and  for  r = 3 
with  sin“l  ^ _ 1+2°,  i.e.,  a = 59° 


o 

9 

f2 

f - f 

il  i2 

f3 

f8 

f9 

0 

0 

0 

0 

0 

0 

0 

10 

.006 

.000 

.006 

.008 

.003 

.111 

20 

.029 

.003 

.026 

.030 

.014 

.222 

30 

.073 

.008 

.065 

.063 

.031 

.331 

1+0 

.11+7 

.016 

.131 

.106 

.0% 

.440 

50 

.258 

.027 

.231 

.156 

.089 

.505 

60 

.1+21 

.01+1 

.380 

.212 

.132 

.649 

70 

.669 

.055 

.614 

.272 

.186 

.745 

80 

1.109 

.068 

1.041 

.334 

.255 

.833 

83 

1.336 

.071 

1.265 

.353 

.278 

.864 

86 

1.693 

.073 

1,620 

.372 

.304 

.878 

89 

2.575 

.075 

2.500 

.391 

.313 

.904 

90 

00 

.075 

00 

.398 

.3^1 

.910 
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TABLE  la 


Values 

of  fc 

>>  r6  and  f7 

for  r = 3 

with 

sin~l  k* 

» 49°,  i.e. 

cc=  59°. 

9 

f4 

4 

f6 

f7 

f -f 
*6  7 

0 

0 

1,000 

-CD 

.419 

-00 

1 

1.015 

-2.537 

.370 

-2.907 

2 

1.031 

-2.072 

.336 

-2,408 

3 

1.047 

-1.793 

.306 

-2.099 

4 

I.O63 

-1.594 

.281 

-1.675 

6 

1.096 

-1.309 

.237 

-1.546 

8 

1.135 

-1.101 

.201 

-1.302 

10 

.011 

1.164 

- .950 

.170 

-1.120 

15 

.023 

1.253 

- .669 

• 113 

- .782 

20 

.037 

1.346 

- .4?8 

.068 

- .54  6 

25 

.053 

1.443 

- .337 

.038 

- .3 75 

30 

.069 

1,544 

- .227 

.020 

- .247 

35 

.084 

1.647 

- .138 

.008 

- .146 

40 

.098 

1,753 

- .064 

.002 

- .066 

45 

.111 

1.862 

0 

0 

0 

50 

.122 

1.974 

.057 

.002 

.055 

55 

.131 

2.088 

.109 

.006 

.103 

60 

.139 

2,205 

.159 

.014 

.145 

65 

.146 

2.326 

.209 

.023 

.186 

70 

.151 

2.451 

.263 

.036 

.227 

75 

.155 

2.580 

.325 

.051 

.274 

80 

.157 

2.714 

.407 

.069 

.338 

83 

.158 

2.797 

.477 

.083 

.394 

86 

.159 

2.882 

.588 

.098 

.400 

89 

.159 

2.970 

.868 

.120 

.748 

90 

.159 

3.000 

00 

.130 

00 
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TABLE  2 


Values 

of  flt  f2 

» f3>  f8  and 

for  r 

= 10 

with 

sin”1  k = 

13°,  i.e., 

a = 72.3° 

<p° 

f 

f 

f 

f 

f 

f 

1 

2 

1 2 

3 

8 

9 

0 

0 

0 

0 

0 

0 

0 

10 

.004 

.000 

.004 

.009 

.006 

.111 

20 

.022 

.003 

.019 

.035 

.027 

.220 

30 

.061 

.010 

.051 

.076 

.057 

.327 

40 

.130 

.023 

.107 

.129 

.103 

.429 

50 

.239 

.041 

.198 

.194 

.158 

.524 

60 

.404 

.064 

.340 

.268 

.230 

.610 

70 

.6  57 

.089 

.568 

.349 

.311 

.685 

80 

1.104 

.114 

.990 

.435 

.406 

.745 

83 

1.334 

,120 

1.214 

.461 

.436 

.759 

86 

1.692 

.124 

1.568 

.487 

.467 

.772 

89 

2.575 

.127 

2.448 

.514 

.498 

.784 

90 

00 

.127 

00 

.523 

.508 

.788 
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TABLE  2a 


Values  of 

V f5>  f6 

and  f^  for  r 

= 10, 

with  sin" 

-1  k'  = 77°, 

i.e.  a = 72, 

.3° 

<p 

f4 

f5 

f6 

f7 

f -f 
7 6 

0 

0 

1.000 

-CO 

1.049 

-00 

1 

.003 

1.063 

-2.424 

.851 

-3.275 

2 

.007 

1.131 

-1.889 

.720 

-2.609 

3 

.015 

1.201 

-1.562 

.607 

-2.169 

4 

.024 

1.275 

-1.329 

.531 

-1.860 

6 

.047 

1.430 

-I.003 

.397 

-1.400 

8 

.073 

1.594 

- .784 

.302 

-1.086 

10 

.099 

1.766 

- .626 

.232 

- .858 

15 

.158 

2.211 

- .379 

.121 

- .500 

20 

.203 

2.667 

- .241 

.063 

- .304 

25 

.238 

3.127 

- .155 

.032 

- .187 

30 

.263 

3.590 

- .097 

.015 

- .112 

35 

.282 

4.054 

- .056 

.006 

- .062 

4o 

.297 

4.523 

- .025 

.001 

- ,026 

45 

.309 

4.996 

0 

0 

0 

50 

.318 

5.478 

,020 

.001 

.019 

55 

.325 

5.970 

.038 

.004 

.034 

60 

.330 

6.476 

.054 

.007 

.047 

6 5 

.33*+ 

6.997 

.069 

.012 

.057 

70 

.337 

7.541 

.085 

.018 

.067 

75 

.339 

8.106 

.103 

.025 

.078 

80 

.340 

8.702 

.12  7 

.033 

.094 

83 

• 34i 

9.074 

.147 

.039 

.108 

86 

.34! 

9.460 

.179 

.045 

.134 

89 

.341 

9.863 

.261 

.054 

.207 

90 

.341 

10.000 

00 

.059 

CO 
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FIG.  4,  S-PLANE 
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CONTRACTION  RATIO  r 

FIG.  5,  COMPARISON  OF  MAXIMUM  WALL  ANGLE 
FOR  DIFFERENT  METHODS 

Present  method  with  <t*  cr„ 

Whitehead,  Wu  and  Waters 

Lighthill 

Cheers 

Libby  and  Reiss 
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CONTRACTION  RATIO  r 


FIG.  5a,  7 AGAINST  r FOR  VARIOUS 

VALUES  OF  a 


FIG.  6,  3‘1  CONTRACTING  DUCT  WITH 


FIG.  7b,  VELOCITY  DISTRIBUTION  ALONG 

DUCT  WITH  a * 72.3' 


